Riemannian geometry of irrotational vortex acoustics 
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We consider acoustic propagation in an irrotational vortex, using the technical machinery of 
differential geometry to investigate the "acoustic geometry" that is probed by the sound waves. The 
acoustic space-time curvature of a constant circulation hydrodynamical vortex leads to deflection 
of phonons at appreciable distances from the vortex core. The scattering angle for phonon rays is 
shown to be quadratic in the small quantity T/{2ncb), where F is the vortex circulation, c the speed 
of sound, and b the impact parameter. 
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Introduction: The last few years have seen considerable 
interest in the development of analog models of and for 
general relativity ||, ^. These analog models provide 
a two-way street: sometimes they illmxiinate aspects of 
general relativity and sometimes the machinery of differ- 
ential geometry can be used to illuminate aspects of the 
analog model. In this article, we use mathematical meth- 
ods developed in the framework of differential geometry 
to study acoustic propagation in an irrotational vortex. 
We analyze the "acoustic" space-time geometry created 
by the vortex in terms of its metric tensor, Killing vec- 
tors, and geodesies. The medium is assumed to be "al- 
most incompressible" , by which we mean that we take 
the background density and the speed of sound relative 
to the medium to be constant, and focus on the effects 
due to motion of the medium. If the flow and its per- 
turbations are irrotational there is a rigorous theorem to 
the effect that sound can be described by a curved-space 
scalar wave equation using an effective acoustic geometry, 
this theorem being derived by using the Euler and conti- 
nuity equations of classical hydrodynamics ||l|, |^. (Even 
if distributed vorticity is present, which is not the case 
considered here, then in the eikonal approximation there 
are rigorous theorems to the effect that sound ray propa- 
gation can be described by curved-space geodesies in the 
same "effective acoustic geometry" Q.) 

Using the acoustic space-time approach, we shall show 
that the quasi-classical scattering process of phonons by 
the vortex leads to a scattering angle quadratic in the 
small quantity r/(27rc6), where T is the vortex circula- 
tion, c the speed of sound, and b the impact parameter of 
the phonon trajectory relative to the vortex. Due to the 
fact that the lowest order scattering angle is quadratic 
in r/(27rc6) and not linear, there is (up to third order) 
no net modification of the transverse force exerted by a 
phonon beam directed towards the vortex. The effect of 



the acoustic space-time curvature is that the vortex acts 
as a converging lens on quasiclassically moving phonons. 

Riemannian geometry of vortex acoustics: For an ir- 
rotational vortex, the velocity profile is given in terms of 
the circulation T = ^ u • dx = 27r 7 by 
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If we adopt cylindrical coordinates r, 9, and z, so that 
the spatial metric is 
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the covariant and contravariant components of the veloc- 
ity are ve — "i] — 7/r^. Now go to a four-dimensional 
description in terms of the coordinates t, r, 9, and z. 
Then, in terms of the speed of sound c, the (3-1-1)- 
dimensional acoustic geometry of the vortex is described 
by the space-time metric [|l|, ^ 
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The scalar velocity potential of phonons lives in the 
curved space-time world given by (?™*°''. This form of 
the metric is often referred to as the Painleve-GuUstrand 
form 1^, 1^ . It is important to note that this vortex met- 
ric is not identical to the metric of a massless spinning 
cosmic string, a solution of the Einstein equations with 
cylindrical symmetry. That metric corresponds to 
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The two metrics (||) and (||), vortex and cosmic string, 
agree asymptotically at large r, 
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Here we have defined the "core radius" 
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This denotes the radius at which the flow goes super- 
sonic. The assumption of incompressibility as well as the 
profile (|l|) certainly break down at this radius, but the 
actual dimensions of the vortex core (for a superfluid vor- 
tex being set by the coherence length ^o) are often larger 
than Tc. The two metrics are markedly different at in- 
termediate values of r. In particular we shall see below 
that they differ significantly well outside the vortex core 
radius Vc- The metric of the spinning cosmic string is 
everywhere flat (the space-time curvature is identically 
zero). In contrast, even at intermediate distances, the 
acoustic geometry of the vortex is not flat (the space-time 
Riemann tensor is not zero), and there are significant ef- 
fects on the propagation of null geodesies (sound rays) 
at values of r well outside the vortex core. That there is 
a region of vanishing classical deflection outside a vortex 
or a magnetic flux tube, like conventionally assumed in 
the standard formulation of the Aharonov-Bohm prob- 
lem (see, e.g., 0|), is thus only asymptotically true 
for a hydrodynamical vortex, even if the generalized vor- 
ticity (hydrodynamical vorticity and/or magnetic flux) 
vanishes everywhere outside the vortex core. It is hence 
not a priori clear that computations of the lordanskii 
force based on the spinning string metric (the analog- 
gravitational Aharonov-Bohm effect (|, ^), give the full 
force exerted by a phonon beam on an actual hydrody- 
namical vortex. 

It is straightforward to compute the Ricci curvature 
scalar corresponding to the metric (H), 
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The curvature of the space-time experienced by the 
quasiparticles is thus significant at distances which can be 
well outside the core domain. That the flow be considered 
well outside the core domain is required by the relation 
(^, which states that at a distance rc from the rota- 
tion axis the velocity around the core equals the speed of 
sound, so that compressibility is no longer negligible. 

A simple criterion for the curvature of the effective 
space-time to be significant for the motion of the phonon 
around the vortex at a given distance r is how the phonon 
wavevector magnitude k compares with the inverse space- 
time curvature radius at that distance. The wavevector 
magnitude is then to be compared with 
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Well outside the core, k^rc ^ 1- The phonon "sees" 
the curvature of the space-time if k exceeds kc{r), and 
behaves as a particle moving on a geodesic in that space- 
time. Vice versa, the topological structure of the acoustic 
space-time related to the lordanskii force (the Aharonov- 
Bohm effect in the space-time of the spinning string) 
dominates the behavior of the phonons if k kc(r). 

For a specific physical example of an irrotational vortex 
geometry, in the dilute limit of a Bosc-condensed atomic 
vapor (BEC) we can relate the circulation F, the 
coherence length and the speed of sound by using 
the relations = h'^ / {"^mgn) and c — y/gn/m. They 
combine to give the relation 

2n rc c = 2TrV2c Co = \T\^27r\-/\ (BEC) (9) 

for a singly quantized vortex with F = h/m, where g 
is the strength of interaction related to the s-wave scat- 
tering length a hy g = inh^a/m. The "acoustic" core 
size is thus in this example of the same order as the 
actual (quantum-mechanical) core size of varying den- 
sity. The hydrodynamic circulation, calculated with the 
speed of sound along a core circumference with radius 
Tc = v^^o, equals the quantum of circulation in the di- 
lute gas limit. This relation qualitatively also holds in 
the dense, strongly correlated superfluid helium II (su- 
perfluid ""He), where is of order the atomic size. The 
relation (||) yields for the curvature scalar in (|^) 
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Curvature of space-time implies that particle world- 
lines on geodesies in that space-time deviate from be- 
ing initially parallel after some proper distance of travel 
along the geodesic. From the general space-time interval 
of the Painleve — GuUstrand metric in the form B 



ds^ = -dt^ + S,j{dx' - v'dt){dx^ - v^dt) , 



(11) 



it is apparent that a constant time slice of the effective 
space-time of quasiparticles is just ordinary Euclidean 
space. Hence spatial distances measured on constant 
time slices in the effective space-time are identical to dis- 
tances in the Newtonian lab world, and a real force is 
acting upon the phonon. It is, however, to be stressed 
that the actual motion of the phonon in the lab world 
is described correctly by the motion in the full effective 
space-time, that is, the phonon is noi just a (Lagrangian) 
particle dragged along by the flow, if that flow is inho- 
mogeneous. 

Phonon motion in acoustic geometry: To explicitly see 
how the vortex flow affects acoustic propagation, and 
thereby get a handle on how acoustic influences can af- 
fect the vortex, we use the eikonal approximation and 
consider phonons instead of sound waves. Phonons then 
follow null geodesies in the acoustic geometry Asso- 
ciated with space-time symmetries are so-called Killing 
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vectors along which the metric is invariant. For both 
geometries (vortex and spinning string) there are three 
such Killing vectors, corresponding to translations in the 
t, 0, and z directions: 



(1,0,0,0); 
and 
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(0,0,1,0); 
: (0,0,0,1), 



(12) 



This leads to three conserved quantities along each 
geodesic [llli 
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Here A is some arbitrary affine parameter for the 
geodesic |ll|. We are interested in null geodesies which 
represent the paths of sound rays in the acoustic geom- 
etry. From these three conservation laws we see that for 
the vortex geometry 
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We are furthermore particularly interested in sound rays 
that come in from infinity, so without loss of generality 
it is possible to re-scale A to choose fci — , and to then 
define fc2 = fc2 and — —k^. Using the new afhne 
parameter, a brief computation yields 



By elimination between the first two equations 
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To now calculate ArjAt we use the fact that the sound 
rays are null curves so that 
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Therefore 
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leading to, by substituting (19) and (|18| 
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The three equations ([T|) ( [T^ , and (|2|) completely spec- 
ify the path of the sound ray in terms of the time pa- 
rameter t and the two nontrivial constants of the motion 
and fcs. These have the physical interpretation that 
^3 = ^00 = /3c < c, while in terms of the impact parame- 
ter h 



j-cb^/l-l3^ 
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The radial motion is more usefully recast as an "energy 
equation" 



1 fArV . 1 2 



with the "potential" 
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The form of this potential is sketched in figure 1 (assum- 
ing ^2 7 > and setting f3 — 0). Note that at large 
distances it has the standard centrifugal barrier propor- 
tional to 1/r^, while at short distances (however, already 
inside the "acoustic" core) it falls quadratically to zero. 




FIG. 1: The effective potential of phonon motion in the 
vortex space-time, with b — —Tc and = v'^ = Tc = c = 1. 

The general deflection angle of the ray as a function of 
/c2 and fcs is obtained by integrating 



r ?dr, 
Ar 
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where the turning point rtum is determined by solving 
dr / dO = 0. For simplicity, we set /3 = fcs = in 
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the following. To lowest order in the small parameter 
\-f\/cb = Tc/b 
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For b — > oo, at fixed 7 and c (fixed core radius rc), the 
deflection then evaluates from integrating 
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along the ray, according to equation (p^). The result is 
to quadratic order in Vc/b 



A0 = 7rsgn(6) 
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The limit of zero classical force corresponds to the limit 
of vanishing core size Vc (c — > cxd and/or 7 — > 0), as may 
be inferred from the vanishing of the curvature scalar (|^). 
This is directly seen in the vanishing, to quadratic order 
in Tc, of the classical scattering of the ray away from a 
straight line caused by the acoustic space-time curvature. 
The action of the vortex on the phonons is schematically 
depicted in figure S. 




FIG. 2: Deflection of phonons by the vortex, which acts as 
a converging lens. 



length is given by 
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In superfluid *He, where — 0.6 A, the roton minimum 
occurs at Ajoton — SA. The phonon wavelength has to 
be about three times the roton wavelength, for an ap- 
proximately linear phonon dispersion to apply, so that 
A 15 Tc, which results in /min — 700 rc. Accordingly, 
the focal length is bounded by / > 40 nm. 

In dilute Bose-Einstein condensates, due to the 
Bogoliubov-type spectrum of these systems, which does 
not display a roton minimum, the phonon dispersion is 
to a good approximation linear up to A ~ 27r^o = V^t^^c, 
hence /min — 20 Tc. With rc — 0.3 ^m, this results in 
f ^ 6 /im for Bose-Einstein condensates. The latter es- 
timate indicates that vortical focusing effects are poten- 
tially within the realm of experimental feasibility. 

Discussion: An incoming phonon of finite momentum 
k passing a singular vortex is deflected by a classical force 
acting upon it. This force is equivalent to an acoustic 
space-time curvature induced in the vicinity of the vortex 
by that flow. From (p9|), it follows that the vortex acts as 
a converging lens. The fact that around a hydrodynam- 
ical vortex there is a classical force field at appreciable 
distances outside the core, entails that calculations based 
on the assumption that there is no force in the vorticity 
free region outside the core are potentially misleading. 
In particular, there are finite distance effects associated 
with a hydrodynamic vortex over and above the familiar 
Aharonov-Bohm effect. 
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From the deflection angle (P9|), one can easily obtain 
the focal length as a function of impact parameter. In- 
deed two initially parallel phonon beams that pass by 
opposite sides of the vortex with impact parameter b will 
intersect at a distance i = 2f beyond the vortex, where 
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(30) 



The fact that the focal length depends on impact parame- 
ter shows that the vortex exhibits what is known in optics 
as "spherical abberation" , which in the present context 
is more correctly termed "cylindrical abberation" . 

In an experiment, one should take into account that 
the transverse size of a phonon beam is limited by the 
wavelength A, which has to be less than or equal to the 
impact parameter b for the quasiclassical phonon scat- 
tering picture to make sense. Thus the minimum focal 
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